Effective mass analysis of Bose-Einstein condensates in optical lattices — 

stabilization and levitation 
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We investigate the time evolution of a Bose-Einstein condensate in a periodic optical potential. 
Using an effective mass formalism, we study the equation of motion for the envelope function 
modulating the Bloch states of the lattice potential. In particular, we show how the negative 
effective mass affects the dynamics of the condensate. 
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I. INTRODUCTION 

The dynamics of classical and quantum systems in pe- 
riodic potentials is a central paradigm of physics, finding 
applications from condensed matter physics to optics, 
and from nonlinear dynamics to atomic and to plasma 
physics. In recent years, the experimental and theoretical 
study of quantum-degenerate atomic systems in periodic 
potentials has opened up new avenues of investigation. In 
particular, Bose-Einstein condensates are a macroscopic 
quantum system that is amenable to exquisite experi- 
mental control. As a result, many phenomena studied in 
solid state systems can be re-examined in a more direct 
and dramatic fashion. Even more importantly perhaps, 
it is now possible to realize experimentally model sys- 
tems that had previously been the object of considerable 
theoretical studies, but were all but impossible to test ex- 
perimentally. One particularly beautiful example is the 
experimental realization of the Hubbard model leading 
to the demonstration of the superfluid to Mott insulator 
transition in a Bose condensate of 87 Rb atoms . 

The first experiments involving the dynamics of Bose- 
Einstein condensates (BEC) in periodic potentials were 
carried out by Anderson and Kasevich, who used this 
approach to demonstrate a mode-locked atom laser Q 
and observe atomic Josephson oscillations ^ |). The 
list of phenomena that were subsequently experimentally 
demonstrated and/or theoretically investigated includes 
the generation of atomic number squeezing Q], the ob- 
servation of the superfluid-Mott insulator phase transi- 
tion, the generation of discrete || or gap solitons Q , the 
prediction of modulational instabilities M and superfluid 
flow H , the observation of Bloch oscillations j| , the anal- 
ysis and observation of coherent acceleration J10[ , studies 
of magnetism etc. 

It is well known that a particle confined to an infi- 
nite periodic potential and acted upon by an external 
force behaves as if possessing an effective mass that can 
be substantially different from its true mass, and may 
even take negative values. In particular, it is this prop- 
erty that is at the core of proposals to generate bright 
matter-wave solitons in BECs with repulsive interactions 
P . The purpose of the present paper is to extend these 



studies by analyzing the temporal evolution of a conden- 
sate in the periodic potential provided by an optical lat- 
tice. The initial state of the condensate is chosen to be a 
(approximate) Bloch state modulated by a slow-varying 
Gaussian envelope, and we compare two situations where 
the Bloch state is either associated with a positive or a 
negative effective mass. We show that in situations where 
the dynamics of the condensate is well approximated by 
a negative effective mass, the periodic potential can lead 
to the stabilization of an otherwise unstable condensate. 
We also demonstrate theoretically the levitation of con- 
densates of negative effective masses. These studies fur- 
ther allow us to determine the impact of self-interactions 
and of finite condensate widths on the general usefulness 
of the effective mass concept. 

The paper is organized as follows. Section II gives a 
brief review of the linear problem of a particle inside a 
periodic sinusoidal potential and introduces the concept 
of effective mass. Section III applies and extends these 
ideas to the dynamics of a BEC in an optical lattice. We 
discuss the equation of motion for the slowly varying con- 
densate envelope function, from which we can gain useful 
physical insights into the dynamics of the system. A full 
numerical solution is presented in Sec. IV, which demon- 
strates in particular the stabilization and the levitation 
of a condensate of negative effective mass, in agreement 
with the prediction of Sec. III. Finally, Section V presents 
concluding remarks on the usefulness of the effective mass 
concept and an outlook. 



II. LINEAR PROPERTIES OF INFINITE 
PERIODIC POTENTIAL 

In this section, we briefly review important aspects of 
the linear problem of a particle of mass m inside a one- 
dimensional infinite periodic potential of the form 



V(x) = V cos 2 (fc x), 



(1) 



with corresponding time-independent Schrodinger equa- 
tion 



- — — + V a cos (M) 



4>{x) = E4>{x). (2) 
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We proceed by introducing the dimensionless quantities 

2mE . o 2m V 



kox; b = 



h 2 k 2 ' 



h A = 



in terms of which Eq. (||) becomes 



d9 2 



h 2 cos 2 i 



h 2 k 2 ' 



0. 



(3) 



(4) 



Eq. (Q) has the form of Mathieu's equation, whose so- 
lutions are well known |l2] |. As a warm-up we sketch 
out the main features of its solutions which, according to 
Bloch's theorem, must be of the form 



.(*) 



isB 



where s is real and arbitrary (in standard textbook lan- 
guage, s = k/ko where hk is the quasi- momentum). In 
our scaled units, the periodic functions F n {9) have a pe- 
riod of 7r. The energy spectrum associated with the pe- 
riodic potential V(6) exhibits a band structure familiar 
from solid state physics. Each value of s gives a discrete 
spectrum whose structure is periodic with respect to s. 
This property allows one to restrict the discussion to the 
first Brillouin zone, — 1 < s < 1. As s is increased within 
that zone, the energy levels trace out curves that are 
restricted to a small energy band. We restrict our dis- 
cussion to the first two bands, as depicted in Fig. 0(a) . 
Note that although the probability distribution \<j)\ is a 
periodic function of 9, <j) itself is not periodic unless s 
takes integer values — e.g., at the center or the edges 
of the Brillouin zone. For a detailed discussion of the 
solutions to Eq. (|J), see e.g. Ref. Q. 

For s = or 1, the eigenstates of the Mathieu's equa- 
tion are periodic and take the form of standing waves. In 
particular, at the edge s = 1 of the first Brillouin zone, 
the wave functions 4> n , s for the first two energy bands, 
n = 1,2, can be expressed as the Fourier series 



!>l, s =l(0) = J2 A2 «+! Sm ( 2 " + 
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(•5) 



(G) 



71=0 



where 
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, n > 1, 



. (a) 



a 2 





1.0 



FIG. 1: (a) Dimensionless energy a = 6 — h 2 /2 and (b) scaled 
reciprocal effective mass m/m* for the first two bands as a 
function of the normalized quasi-momentum s. In (b), the 
solid and dashed lines refer to the first and second band, re- 
spectively. For this calculation, h 2 — 3.2. 



The effective mass, to*, which characterizes the re- 
sponse of the particle to external perturbations, is defined 
as 



2m 



d 2 E/dk 2 d 2 b/ds 2 



(7) 



Fig. |l|(b) shows the variation of the effective mass of the 
first two bands with the quasi-momentum. One observes 
that it is near infinite (i.e., 1/m* = 0) for certain values 
of the quasi-momentum s, and even becomes negative in 
certain regions of the zone. In more than one dimension, 
to* acquires a tensorial character, with elements given by 



h 2 



"fiv 



d 2 E/{dk^dkv 



In the next section, we apply the effective mass formal- 
ism to the case of a self-interacting Schrodinger field de- 
scribed by the Gross-Pitaevskii equation. Specifically, we 
show that a BEC prepared in a Bloch state with negative 
effective mass behaves as if the signs of its self- interaction 
and of the external confining potential had been reversed. 



k n = (46- 2h 2 -An 2 )/h 2 



III. EFFECTIVE MASS EQUATION 



We consider a BEC in a sinusoidal optical lattice 
G 2n +i = A2n+i/A 2n -i, or B 2 n+i/B 2 n-i, potential. Near zero temperature, the system is de- 

scribed to an excellent degree of approximation by a time- 
and A3 = fei + 1, B3 = fei — 1. dependent Gross-Pitaevskii equation for the normalized 
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condensate wave function ^>(r, t), generalized to include the effects of three-body collisions, 
I 



ift 



dr />(T,t) 
dt 



_ v 2 + V(r) + U(r) + Ng\ij\ 2 - iN 2 K 3 \^\ 4 
2m 



(8) 



r 



Here, V(r) is the periodic lattice potential, U(r) is an 
additional external potential that is taken to be slowly 
varying on the scale of the lattice period (e.g., a confin- 
ing potential), TV is the total number of atoms, g is the 
strength of the nonlinear inter-atomic interaction, and 
K 3 gives the rate of three-body recombination loss E3] . 

We proceed by expanding the condensate wave func- 
tion ijj(r, t) on the complete set of Bloch functions n k(r) 
that are the stationary solutions of the associated linear 
Schrodinger equation 



2m 



V(t) 



</>«k(r) = -E„k</>nk(r). 



(9) 



The subscripts n and k represent the band index and 
quasi-momentum, respectively. The Bloch functions sat- 
isfy the orthonormality condition: 



dr Cik( r ) </w (r) 



kk' 



where the integral is taken over a single period of the 
lattice of volume Q. Our use of the resulting expansion, 



tp(r, t) = A„ k (r, t)<j) nk (r)( 



-iE nk t/h 



(10) 



is motivated by the the fact that the Bloch functions 
can capture the rapid oscillations of the condensate wave 



J 



function, while the slowly varying envelope functions A n k 
will describe the slow center-of-mass motion of the con- 
densate. 

The equations of motion for the amplitudes A n ^ are 
found by inserting the expansion ( |Io| ) into the general- 
ized Gross-Pitaevskii equation (^). The resulting equa- 
tions are however of little use if a very large number of 
A n k is required to accurately describe the dynamics of 
the condensate. Here, we consider a more specialized sit- 
uation where the matter-wave field is characterized by 
a central wave vector ko corresponding to the mean ve- 
locity of the condensate (this would correspond to the 
carrier wave in conventional optics). Hence, we expand 
■0(r, t) in a way reminiscent of the slowly varying enve- 
lope approximation of quantum optics as 



V>(M) = fn{r,t)(j) nko (r)e- 



iE nk() t/h 



(11) 



Inserting Eq. (|llj) into Eq. 
mass jL4| or multiple scales 



and applying the effective 
15[ methods yields a set of 
equations governing the time evolution of the condensate 
envelope functions f n (r,t) 



ih 



dfn 

dt 



if- 



2m* dx^dxv 



+ U(r) + Ng'\f n \ 2 -iN 2 K' 3 \f n 



(12) 



r 



The derivation of Eq. (jT|), the so-called effective mass 
equation, is given in Appendix |A|. Here x^ and Xu are 
cartesian components of r and we use an implicit sum- 
mation over repeated indices. The velocity 



*>nko|p|0nk o ), 



(13) 



is the drift velocity of the n-th band contribution to the 
condensate envelope. As shown in Appendix |A|, it is equal 
to the gradient of the energy _E„k with respect to the 
quasi-momentum k, evaluated at kp. We note that v„ 



vanishes at the extreme points of the band, in particular 
at the band edges. The coefficients g' and K' give the 
renormalized nonlinear interaction strength and three- 
body loss rate, respectively. They are given by 



and 



Q 



g / dr |0 nko (r)| 4 , 

J cell 

K 3 [ dr|0„ ko (r)| f 

./cell 



(14) 



(15) 
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The lattice potential V(r) does not appear in these ex- 
pressions, its effects being already incorporated into the 
effective mass to* which appears in the first term on the 
right-hand side of Eq. ( jl2| ) as well as in the drift term 



J 



on the left-hand side of that equation. 

The effects of a negative mass, to* < 0, can be imme- 
diately inferred from the complex conjugate of Eq. (|l2|), 



2 1 to* I dx^dx l 



-U(v)-Ng>\r n \ 2 -iN 2 K'z\f n \ 



f* 



(16) 



This equation shows that as compared to the case of a 
positive effective mass, the condensate density envelope 
profile |/n(t)| 2 evolves now under the influence of an in- 
verted external potential —U(r) and an inverted non- 
linearity —g'. As expected, though, the three-body re- 
combination rate does not change sign and still repre- 
sents a loss term. The change in sign of g 1 and U leads 
to a number of consequences. In particular, past work 
has shown how this property can be exploited to launch 
bright matter-wave solitons in condensates with positive 
scattering length || and induce modulational instability 
0. The following section explores further consequences 
of this property. 



IV. EFFECTS OF NEGATIVE EFFECTIVE 
MASS 

A. Inverted nonlinearity 

We first investigate the effect of an inverted nonlin- 
earity. As is well known, this will change the two-body 
interaction from being repulsive to attractive, and vice 
versa. Since condensates with attractive interactions are 
normally unstable against collapse, the sign change of to* 
is clearly expected to have a significant effect. 

To reduce the scale of the computation, we restrict 
ourselves to a two-dimensional system consisting of a 
pancake-shaped condensate whose axial motion is frozen 
to the trap ground state. (We note that condensates 
with reduced dimensions have already been realized in 
the laboratory |fl7| .) Furthermore, to focus on the effects 
of the nonlinearity, we assume that except for the lattice 
potential 

V(x, y) = V [cos 2 (k Q x) + cos 2 (k y)), (17) 

there is no additional external potential in the transverse 
dimensions. 

The evolution of the transverse component ip± of the 
condensate wave function is then given by the generalized 
Gross-Pitaevskii equation 

_dtp± 



where 



ih- 



-— VI + V(x,y) + Ng\^ ± \ 2 - iN K 3 \ip±\ 
Zm 



Vi = d 2 /dx 2 + d 2 /d v 2 



For the initial state of the condensate, we assume a gaus- 
sian transverse wave function of the form 



i>±(x,y,t = 0) =Aff(x,y) c/) nk (x,y), 



(18) 



where f{x,y) = cxp[— (x 2 + y 2 )/w 2 } is a slowly varying 
Gaussian envelope function of width w, 4> n k is a Bloch 
wave function of the linear Hamiltonian, and the con- 
stant Af normalizes ip± to unity. Figs. || and || illustrate 
the dynamics of the condensate for two cases, one cor- 
responding to a negative effective mass and the other to 
a positive one. We choose the Bloch wave function 
associated with a negative effective mass to* as 



0„ k ~ sin(fc a;) sm(k y). 



(19) 



This corresponds approximately to the Bloch state at 
the top of the first energy band, whose full expression 
is given by <^i )S =i of Eq. (||). The approximate form 
19| ) retains only the dominant term in the Fourier series 
. Its overlap with <^i, s =i is larger than 99% for the 
parameters of our simulations. Similarly, for the state of 
positive effective mass to* we use 



4>nk ~ COS(/C X) COs(fc 2/), 



(20) 



which corresponds roughly to the Bloch state 4>2, s =i at 
the bottom of the second band, see Eq. @ . Both states 
@H) and @ are located at the edge of the first Brillouin 
zone, with corresponding drift velocities v g = 0. Hence 
the center of the condensate stays at x = y = during 
the course of its time evolution. 

The two columns in Fig. || show snapshots of the den- 
sity distribution of the condensate for positive (left col- 
umn) and negative (right column) effective masses to*. 
In both cases the nonlinear interaction strength g is neg- 
ative, hence in the absence of the lattice the condensate 
would be unstable and subject to collapse. This how- 
ever can be changed if the effective mass of the BEC 
becomes negative, as demonstrated in the right column. 
In that case, the condensate behaves as if experiencing 
a repulsive two-body interaction, its width increasing in 
time. This is in sharp contrast to the situation for a pos- 
itive effective mass, in which case the condensate width 





FIG. 3: Temporal evolution of the number of atoms in the 
condensate. The solid (dashed) line is for the positive (nega- 
tive) m* state. Same parameters as in Fig. 0. 





FIG. 2: Temporal evolution of the density distribution of a 
BEC in optical lattice. Lighter shades of grey correspond to 
higher density. The left (right) panel shows the evolution of a 
positive (negative) m* state. The dimensionless time is t~0, 
1, 2, and 3 from the top to bottom. In the plots, the ranges for 
x and y are from —30 to 30. The condensate parameters are 
N = 50000, g = -0.004, V = 10, fc = 3, K 3 = 3.9 x 10~ 8 
and the width of the initial envelope function w = 8. The 
units for time, length and energy are 1/u, y H/Jnuj), and 
hu, respectively, where ui is the trap frequency in the axial 
direction. 



quickly implodes jL8|. This implosion is characterized 
by a complex dynamical behavior, such as the forma- 
tion of atomic "bursts," observed at later times. (We 
remark that while an attractive condensate can expand 
if it possesses a sufficiently large kinetic energy, this is 
not the case of Fig. ||, where the two initial states have 
nearly identical kinetic energies.) The collapse dynamics 
depicted here is rather similar to that of an attractive 
condensate in a harmonic potential, a system which has 
gained much interest recently as a result of the capability 
of tuning the scattering length from positive to negative 
via Feshbach resonance |2(], 21 



We note however 

that the use of the generalized Gross-Pitaevskii equation 
becomes questionable in the regime of extreme collapse. 
This point though it is important, is not central to our 
results. 

We conclude this section by observing that the number 



of condensate atoms decreases in time, due to the pres- 
ence of three-body loss. This is particularly evident for 
the state with positive effective mass, as shown in Fig. ||. 
This is as expected, since the atomic density is higher in 
that case. Finally, we note that as a result of the tensorial 
character of the effective mass, it can readily be positive 
along some dimensions and negative along others. For 
example, in two dimensions a condensate wave function 
of the form 

0„ k ~ sin(fc a;) cos(fc y), 

has an effective mass that is negative along x and positive 
along y-direction. This leads to a situation where the 
condensate will expand along x and contract along y, as 
has been numerically confirmed. 



B. Inverted external potential 

We mentioned that if a condensate of negative effective 
mass is subject to an external potential U(r) that varies 
slowly over the lattice period, then this potential will 
appear inverted to the condensate. For example, under 
the effect of gravity, the center of mass of the condensate 
will climb up, rather than fall down. To demonstrate 
this effect, we solve the one-dimensional Gross-Pitaevskii 
equation 



ih 



dip 

"at 



2m dz 2 



+ V(z) + U(z) + Ng\*p\ 2 



4>, (21) 



where V(z) = Vo cos 2 (koz) is the lattice potential and 

U(z) = Gz, 

represents the gravitational potential, G being a con- 
stant. We neglect the three-body loss term, which is 
not essential in the present discussion. 

Again, we assume that the initial condensate wave 
function is a broad Gaussian envelope function centered 
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FIG. 4: (a) Temporal evolution of the condensate center of 
mass, (b) Solid line: density distribution of the state with 
negative effective mass at the dimensional time t — 7.53 when 
the center of mass of the condensate reaches maximum height. 
The dashed line shows the initial Gaussian envelope function 
for reference. The parameters used are Vo = 15, w = 8, 
ko = 5, G — 0.65, and Ng = 10, in the same dimensionless 
units as in Fig. is used. 



gravitational energy. For lithium atoms in an optical 
lattice formed by laser light of wavelength of 0.5 /xm, 
one finds that z ma x is about 0.4 mm, a rather significant 
distance in the context of integrated atom optics [Bp . 
Figure ^(b) compares the condensate density distribution 
at the maximum height to its initial envelope profile. 

That a maximum height should exist demonstrates the 
limit of the effective mass formalism, since a literal in- 
terpretation of the effective mass equation ( |l2|) would 
incorrectly lead one to conclude that in the case of neg- 
ative effective mass, the condensate should move against 
the gravitational field without bound. The incorrectness 
of this reasoning lies in the fact that Eq. ([l2|) neglects 
all terms coupling different Bloch modes (see derivation 
in Appendix [A|). What happens in reality is that as 
the condensate moves under the combined influence of 
the external potential and of its intrinsic nonlinear po- 
tential, other Bloch modes with different quasi-momenta 
than those initially excited inevitably become populated. 
This in turn changes the effective mass of the system, or 
even makes this band-specific concept meaningless. In 
the present example, the effective mass changes its sign 
to positive at the maximum height as successive energy 
bands become excited. In the mean time, the oscilla- 
tions in the density profile have a contrast less than one 
[see Fig. [|(b)] which shows that the phase of the wave 
function at maximum height becomes spatially continu- 
ous (otherwise, the density would vanish at each phase 
singularities) while the initial wave function possesses a 
series of 7r-phase jumps. 



C. Discussion 



at z = 0, modulated by a sinusoidal function. Taking 
that modulation of the form sin(fcoz) yields a state of 
negative effective mass m*, while cos(fcoz) gives a state 
of positive effective mass. Figure |I](a) illustrates the evo- 
lution of the center of mass of the condensate for these 
two cases. As expected from the previous discussion, we 
see that it initially climbs up for negative m* and falls 
down for positive m* . 

As the center of mass of the condensate of negative 
effective mass climbs up against the gravitational field, 
its gravitational potential energy increases. This increase 
must be offset by a decrease in kinetic energy (assuming 
for the sake of argument that the nonlinear interaction 
energy is negligible). Since the initial condensate pos- 
sesses a finite kinetic energy, roughly given by fi 2 fcp/(2ra), 
conservation of energy therefore sets an upper limit on 
how high its center of mass can climb. This is illustrated 
in Fig. ^(a). A good estimate for the maximum height 
can be obtained from 



where the left hand side is about half of the initial kinetic 
energy and the right hand side is the maximum gain in 



We conclude this section with some general comments 
on the generation and robustness of condensate states of 
negative effective mass. As one can see from Eq. (fTs|), 
the condensate wave function suffers a 7r-phase jump ev- 
ery lattice period. This immediately suggests that these 
states can be prepared by the phase-imprinting method 
that has been previously successfully implemented to 
generate dark solitons as in Ref. (^3). However, due to 
edge effects and inaccuracy in controlling the intensity 
of the phase imprinting pulse, this method may not be 
fully adequate to generate a series of accurate 7r-jumps 
as required by Eq. ( p^ ) . 

A better choice may therefore consist in first prepar- 
ing a condensate with uniform phase in a specific elec- 
tronic state. Taking advantage of the internal degrees 
of freedom of the atoms, one can then drive the conden- 
sate into another internal state with an optical field of 
appropriate topological spatial geometry. The resulting 
center-of-mass wave function of the condensate having 
then the correct phase structure. This method has been 
successfully used in the past to create a series of accurate 
7T- phase jumps in condensates |^4|, |25|, |2r| . 

We conclude this section with a brief discussion of the 
robustness of condensate states of negative effective mass. 
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As we have mentioned, the state (|T|) is only an approxi- 
mation to the true Bloch state (j^) . The overlap between 
these two states decreases as h 2 increases. Furthermore, 
|l/m*| decreases with increasing h 2 , as the system ap- 
proaches the tight-binding regime. On the other hand, 
h 2 cannot be made too small, since the size of the gap 
between the first and second energy band is itself pro- 
portional to h 2 , and a small gap will result in inter-band 
transitions to a state of positive effective mass. The op- 
timal value od h 2 also depends on other system parame- 
ters, and as a result we found that h 2 ~ 2 turned out to 
be optimal for our purpose. 



V. CONCLUSION 

In summary, we have studied the dynamics of an 
atomic condensate in a periodic optical lattice using an 
effective mass method. We have obtained an effective 
equation of motion governing the time evolution of the 
envelope of the condensate wave function in which the 
periodic external potential appears in the form of an ef- 
fective mass, which can be either positive or negative, 
and a global drift. Numerical calculation confirmed that 
this envelope function approach provides us with useful 
qualitative insight into the condensate dynamics. Our 
study focused on negative effective masses, in which case 
the condensate behaves as if it subject to an inverted 
nonlinearity and an inverted external potential. 

We emphasize that caution must be used when apply- 
ing the effective mass concept, which neglects the poten- 
tially important coupling between different Bloch modes. 
This is particularly so for a system with negative effec- 
tive ma ss mo ving in an external potential, as illustrated 
in Sec. IV B. An ordinary particle inside an external 



potential will move towards the potential minimum, in 
the mean time gaining kinetic energy. A particle with 
negative effective mass, however, will move towards the 
potential maximum and in doing so, it looses kinetic en- 
ergy. This process cannot go on forever since there is 
only a finite amount of kinetic energy for the system to 
loose. In the process, the effective mass will eventually 
change from negative to positive. 

Our study also demonstrates that, as an alternative 
to the established Feshbach resonance method, changing 
the effective mass of a Bose-Einstein condensate provides 
us with an additional way to control the nonlinear atom- 
atom interaction in a condensate. This approach should 
be particularly useful for systems with no Feshbach res- 
onances at convenient magnetic field strengths, or when 
the presence of external magnetic fields is undesirable. 
Controlling atom-atom interactions via their effective 
mass also provides the possibility to induce anisotropic 
nonlinear interactions, due to the tensorial character of 
the effective mass in higher dimensions. 

Our study here focused on the dynamics of the sys- 
tem. In the future, it will also be interesting to study 
its static properties. For doing so, one needs to provide 



an additional confining potential. We must keep in mind 
that a confining potential for a negative effective mass 
state is an anti-trap, instead of a trap. This suggests the 
possibility of trapping strong-field-seeking states using 
the magnetic traps implemented in current BEC experi- 
ments. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE MASS EQUATION 

In this Appendix, we give a derivation of the effective 
mass equation (12). For the sake of simplicity, we restrict 
ourselves to a one-dimensional system and neglect the 
nonlinear terms, which can be added straightforwardly. 

We are concerned with the solution of the Schrodinger 
equation 



dj>{z,t) 

at 



= [H + U(z)}TP(z,t), 



(Al) 



where H = p 2 /2m+V(z) is the Hamiltonian for periodic 
potential V and U is some additional potential that varies 
on a spatial scale much larger than the lattice period. 



1. The k ■ p perturbation method 

Associated with each quasi- momentum hk is a set 
of Bloch functions ^nk with energy E n ^ as defined in 
Eq. (^) . In this section, the power expansion around the 
central wave vector kg of the energy band function E n k is 
investigated. In particular, we seek expressions for first 
and second derivatives of E n k with respect to fc. These 
expressions will become useful in the derivation of Eq. 

©• 

The "k • p" perturbation approach is a straightfor- 
ward method to relate the energies and wave functions at 
nearby points in the quasi-momentum space. We proceed 
by writing the Bloch function </> ra k as 



4>nk(z) 



e lkz u nk (z), 



where u n k is the cell periodic function. From Eq. 
immediately have 



p 2 | | h 2 k 2 

2m 2m 



hk ; 

— P 

m 



), we 



(A2) 
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We then perform a power expansion around fco up to 
second order, such that 



k = fcg + 5k, 



1 ,„. ,r> A d 2 E n k 



dk 2 



ink = u nka + {6k) 2J a^^mfco + (<5fc) 2 ^ a;,: '"/„/.„• 



,( 2 ), 



where (■ ■ • )o means that the corresponding derivative is 
eval uate d at k = fco- Next we insert these expansions into 
Eq. (A2), and equate the coefficients of different powers 
of 8k. For the zeroth order, we find 



A 2 
P 

2m 



V(z) 



h 2 k 2 

2m 



hkn 



-P 



t^nfco — E nko Unko 



which is simply the definition of u nko . The first-order 
term yields 



Since 1/m* = (d 2 E nk /dk 2 ) /h 2 , we have 



1 1 1 



\ " 



m* m m 2 *-f E 7nko E nko 



(A4) 



2. Derivation of Eq. (12) 



We now proceed to derive Eq. (|lj, adaptingwith slight 
modifications from the approach of Ref.jl4|]. First, we 
introduce the Fourier expansion of the envelope function 
as 



f n (z,t)= / A nk (t)e^ k - k ^ dk. 



fdE n 



n 2 k 



\ dk J m m 

^ t a rn {Emko ~ E nko )u mka . 



Unkn 



(A3) 



Taking the dot product of Eq. ( |A3| ) with u nko and inte- 
grating over z in the first Brillouin zone, we find 

1 ( dE nk \ hk 1 

"FT, — = 1 Pnn{k ) 

ok J m m 

Kk \P\4>nk ) = V g , 



1 



From Eq. (pd|), we have 



i/>(z,t) = E/ A nk (t) X nk(z)e- lE ^o^dk, (A5) 



where 



Xnk(z) = e l(fe ~ feo) ^„ feo (z) = e lk *u nko (z). 



where we have used the orthonormality condition of the 
Bloch functions and 



After some algebra, one finds the matrix elements of 
Hq and U as 



2tt 

Pnl(k) = —(u nk \p \uik). 

One can see that the drift velocity as defined in Eq. Jl3| ) 
is just the gradient of the band energy at fco- 

Carrying out a similar procedure with u mko (m ^ n), 

we get from (A3) the expression for Om as 



(Xnk\H \Xlq) 



S(q-k) 



E, 



ik„ 



h 2 (q 2 -k 2 ) 
2m 



Snl + 



H(q - fc ) 



in 



Pnl(k ) 



(Xnk\U\xi q ) = S nt U(k - q), 



Pn 



m E mka — E nko 

Analogously, by using the second-order equation, we can 
find the second derivative of E nk as 



where 



U(k - q) = — ( U{z)e l{q - k)z dz. 
2tt J 



d 2 E n ] { 
dk 2 



,2 



m m 



Emko E nko 



Inserting (|A5|) into (Al), we have 



ih 



8A nk (t) _ h 2 (k 2 -k 2 ) A t h(k-k Q ) 



dt 



2m 



nk 



p n i{k )A lk + / dqU(k - q) A nq . 



(A6) 
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The only terms that couple different bands (I ^ n) are 
those involving the momentum matrix elements p n i. As 
long as these inter-band couplings are weak enough, we 
can treat them as weak perturbations. To lowest order, 
we have 



.4 



Ik 



h(k-k ) pi n (k ) 



■ A r , 



m Eik — E nko 
Inserting this expression into Eq. (|A6|), we have 



. dA nk 



h 2 {8kf 

2m* 



Hv g (5k) 



A nk + / dqU(k-q) A nq , 



where 8k = k — kg. Performing an inverse Fourier trans- 
form from A n k(t) back to f n (z,t) and replacing 5k by 
its configuration space representation —id/dz, we finally 
obtain the one-dimensional version of Eq. ( |l2| ) . The non- 
linear terms can be added to this equation in a straight- 
forward way. 
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